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For quantum mechanics on a lattice the position ("particle number") operator and the quasi- 
momentum ("phase") operator obey canonical commutation relations (CCR) only on a dense set 
of the Hilbert space. We compare exact numerical results for a particle in simple potentials on the 
lattice with the expectations, when the CCR are assumed to be strictly obeyed. Only for sufficiently 
smooth eigenfunctions this leads to reasonable results. In the long time limit the use of the CCR 
can lead to a qualitatively wrong dynamics even if the initial state is in the dense set. 



Since the pioneering work on_quantization of the elec- 
tromagnetic field, where Diractl introduced a phase ob- 
servable, supposedly conjugate to the number operator 
Af, there has been a long controversy wether a phase op- 
erator 9 can be constructed which obeys the canonical 
commutation relation (CCR) [Af, 9] — il. AsJiere exist 
excellent reviews which address this questiorou it is not 
necessary to give a detailed history of the various propos- 
als. We want to point out that there are different types 
of "number operators" of importance. In quantum op- 
tics the eigenvalues of the particle number operator run 
from zero to infinity. For the case of fermions with a filled 
Dirac sea one often works with the number operator rela- 
tive to the Dirac sea which has integer eigenvalues which 
run from minus infinity to plus infinityuu. A concep- 
tionally simple case where the number operator is of the 
form Af — X^oo l n ) n ( n l is (apart from the multiplicative 
factor of the lattice spacing a) the position operator x for 
quantum mechanics on a one-dimensional lattice!!. In this 
paper case we treat infinite lattices which corresponds to 
the simpler case of a number operator not bounded from 
belowta, while the quantum optics scenario corresponds 
to a semi-infinite lattice. 

Shortly after Dirac's introduction of the phase opera- 
tor Jordan] gave an argument why the CCR [Af, 9} = il 
cannot hold as an operator identity. This commutation 
relation would imply [sin aAf, 9] — ia cos cxhf. If one 
puts a = 7r this leads to a contradiction as sin7rjV is 
zero because Af has integer eigenvalues, but cos7r7V is 
nonzero. It is only possible to construct a canonical jp. 
(Af, 9) on a dense set of the states of the Hilbert spaceB'E. 
The question we want to address in this paper is if this 
fact is useful for typical applications in physics, where a 
Hamiltonian is given and one wants to solve the time- 
independent and (or) the time-dependent Schrodinger 
equation. An example where in this context the canonical 
commutation relation [Af, 9) = il is generally-assumed is 
in treatments of the Josephson HamiltonianHl. 

As a simple example we study the properties of a par- 
ticle described by the Hamiltonian 



H = - ^2(\m + n)t n (m\ + H.c.) + \m)V„, 



(1) 



where the t n are the hopping matrix elements and the 
V m are the local potentials. For a constant potential the 
(improper) eigenstates are the Bloch states 



(2) 



with k in the first Brillouin zone (—n/a, Tt/a]. The Bloch 
states are normalized as (k\k') = 6(k—k'). The "obvious" 
operator to replace the momentum operator (divided by 
H) on the lattice is the quasi-momentum operator 
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k= \k)k{k\dk. 

J —-K I a 



(3) 



In the position representation 9 = ak is just the infinite 
lattice version of the "phase operator" proposed inde- 
pendently by Garrison and Wongp and GalindoEJ for the 
semi-infinite case 
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If one applies x = aAf to the states in Eq. (4) one sees 
that x9\n) is not an element of the Hilbert space li, while 
9x\n) is. A straightforward formal calculation shows that 
x and k do not obey a CCR 



x , k 
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The correction term can be expressed in tetms of the 
projection operator on the Bloch state |7r/a)c3 
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For states \ip) where the components ip n = (n\ip) obey 
the condition = J2n(~ — tne correction term 

vanishes. These_s.tates form a dense set of states of the 
Hilbert space ijEi. This condition is fulfilled when the 
components in position space are sufficiently smooth. If 
4>n = f{an), where / is a differentiable function, the 
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operator ik for a — > acts like the differential operator 
d/d(an). Using Eq. (4) this can be seen by writing 

(n\ik\tp) = 22 : i^n+j - ipn-j) ■ (7) 

i=i 3 

With the assumption on the function / made above the 
series on the right hand side (rhs) is approximately of 
the form 2(/'(an) — f'(an) + f'(an) — h..) which quickly 
converges to /'(an)..— .This can be best seen using the 
Euler transformation!!] of the sum on the rhs of Eq. (7). 

We next consider junctions of the operator k. For in- 
teger n the exponential functions 

e -mak = ^2 \ m + n) (m| = f n (8) 

m 

are just the discrete translation operators. Using the site 
representation one easily shows the commutation relation 
[x,e~ mah ] = ane~ mah , which also can be "derived" as- 
suming the CCR. The latter incorrect derivation could 
also be performed if n is replaced by an arbitrary real 
number. But then the above result no longer holds, i.e. 
the use of the CCR gives a wrong answer. In the follow- 
ing we also use the operator k . Its matrix elements in 
the position representation are given by 

7r 2 (—\\m-n 
a 2 (m\k 2 \n) = —6 mn + 2 '- (1 - S mn ). (9) 
6 (to — n) z 

Therefore the operator k 2 /2 has the form of the first 
term on the rhs of Eq. (1) with to = — 7r 2 /(6a 2 ) and 
t n = (— l) n+1 / (an) 2 . The corresponding energy eigen- 
values efc = — to — 2^2 n>l t n cos (akn) are just the well 
known Fourier series which corresponds to the periodi- 
cally continued parabola arcs. In the following we mainly 
work with the Hamiltonian H = k /2 + V(x), which is a 
special case of Eq. (1). 

We examine, if assuming the strict validity of the CCR 
[x,k] = il can be useful in the description of a particle 
governed by such a Hamiltonian by comparison with ex- 
act results. We study the eigenvalue problem as well 
as solutions to the time dependent Schrodinger equa- 
tion. Only for those eigenstates \ipE) which approxi- 
mately fulfill (Tr/a\ipE) = Owe can expect the assumption 
[x, k] = il to be useful. 

We begin with the harmonic oscillator on the lattice 
H = fc 2 /2 + ex 2 /2_ The Josephson Hamiltonian in its 
most simple formEll correponds to a quadratic poten- 
tial but with a kinetic energy Tki n with nearest neigh- 
bour hopping (1 - cos ak)/a 2 = [i - (T x + T_i)/2]/a 2 
for a = 1. We present results for both types of the ki- 
netic energy and explicitely mention when T^m = k 2 /2 
is not used . The exact solution can be obtained numer- 
ically by diagonalizing a sufficiently large matrix. The 
expectations from assuming the CCR [x,k] = il are ob- 
vious if one introduces the usual ladder operatorsQ and 



determines the eigenvalue spectrum purely algebraically, 
which yields E n = y/c(n +1/2). If we further assume 
the expectation value of the potential energy to be half 
the energy eigenvalue as in the continuum limit it follows 

(n) 

that the the wave function ij)m approximately extends 
from —(n/^/c) 1 / 2 to (n/^/c) 1 / 2 . The mean separation of 
the nodes of the corresponding wave function in the con- 
tinuum is therefore given by 2/Vnc 1 / 2 . The smoothness 
assumption which is necessary for the CCR to be approx- 
imately valid certainly breaks down, when this separation 
becomes of the order of the lattice constant a. One there- 
fore expects the equidistant energy spectrum of the har- 
monic oscillator in the continuum to be app roxim ately 
valid for quantum numbers n smaller than h/ V a 4 c where 
b is of order unity. The exact numerical results for the 
lattice model are shown in Fig. 1. 
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FIG. 1. Normalized energy eigenvalues of the harmonic os- 
cillator on the lattice as a function of ac 1//4 . Below the 
dashed line 3/(ac 1/ ^ 4 ) 2 explained in the text one approxi- 
mately recovers the eigenvalue spectrum of the continuum 
limit. The dotted line shows two eigenvalues for the kinetic 
energy (1 — cosafc)/a 2 . 

The dashed line shows the criterion just derived for 6 = 3. 
Above this line there are strong deviations from the 
equidistant spectrum, in fact pairs of eigenvalues become 
almost degenerate. A simple explanation for this obser- 
vation will be given later. For comparison we also show 
two eigenvalues for the nearest neighbour kinetic energy 
(1 — cosak)/a 2 which for small enough values of ac 1 / 4 
are expected to form the same equidistant spectrum. As 
shown the results deviate from E n = \fc\n + 1/2) al- 
ready for smaller values of ac 1 / 4 and the almost degen- 
eracy of the eigenvalue pair occurs already for a smaller 
value of ac 1 / 4 compared to the long range hopping. Fig. 
2 shows the behaviour of the absolute value of the overlap 
S n = 53 m (— l)ipm as a function of the quantum number 
n for the even parity states. The values of n where S n 
starts to become "nonzero" agree very well with expec- 
tation from the argument using the number of nodes. 
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In order to understand the spectrum in the upper right 
corner of Fig. 1 it is useful to first treat the case of a linear 
potential described by the Hamiltonian H = k jl — Fx. 
For such a potential the energy eigenvalue spectrum for 
the lattice model is known to be completely different 
from the continuum limit. It consists of an equidistant 
"Wannier-Stark"-laddeili-a, while the spectrum is purely 
continous for the continuum case. This can easily be 
understood using the translation operators T n defined 
in Eq. (8) which obey the simple commutation rela- 
tion [T n , H] = anFT n . After one localized eigenstate 
\Eq) with eigenvalue Eq is shown to existH-3 this commu- 
tation relation implies that the T n \Eo) are eigenstates 
with eigenvalues Eq — anF. The form of the localized 
Wannier-Stark states for the long range hopping of Eq. 
(9) is shown in Fig. 3 (stars). The other wave function 
(circles) is discussed when we return to the discussion of 
the harmonic potential. 
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FIG. 2. Absolute value of the overlap S n = ^2 m {—l)'4>m > 
which determines the deviation from the CCR as a function of 
the quantum number n for the even parity states, and different 
values of c: c = 1 squares, c = 0.1 stars and c = 0.01 circles. 

We next examine if the CCR assumption is of any use 
for the understanding of the time evolution in the linear 
potential. The equations of motion for the Heisenberg 
operators £(t) = e lHt xe~ lHt and k (t) close if the CCR 
are assumed and yield fc(£)cCR = k + Ftl as well as 
x(t)ccR — x + kt + lFt 2 /2. Obviously the result for 
k(t)ccR is incorrect as the operator k defined in Eq. 
(3) is bounded and this property also has to hold for 
k(t) . The result for the position operator using the 
CCR differs from the exact solution, which follows from 
[x, T n ] = anT n and [T„, H] = anFT n as 



i(t) =£ + 
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- l) jF + H.c. 



(10) 



This shows that the exact solution is periodic in tipp. 
The particle undergoes so called "Bloch-oscillations"o 
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FIG. 3. The stars show the amplitudes v2Vm of the Wan- 
nier-Stark eigenstate for a linear potential on a lattice with 
F = 0.4 localized around the site m = —41. The circles show 
the amplitudes ip m of the even eigenstate of the harmonic 
potential with c = 0.01 localized around the same site. 

The two results are compared in Fig. 4 for two Gaus- 
sian wave pakets (ip n (t = 0) ~ exp [— b(n — no) 2 ]) with 
the same location but different widths as initial states. 




FIG. 4. The full curve (6 = 0.2) and the dashed-dotted 
(b — 0.02) curve show the periodic Bloch oscillations in a 
linear potential (F = 0.4). The long dashed curve shows the 
result from assuming the CCR. In order to understand why 
this approximation fails qualitatively when approaching the 
time ir/(aF) the time dependent overlaps IS^t)! (see text) are 
shown as the dotted (b = 0.2) and dashed curves (b — 0.02) . 

If the initial wave paket extends over many lattice sites 
the CCR assumption yields an excellent approximation 
for times less than Tt/(aF). At this time the turning 
point is reached in the exact solution. Approaching this 
time the overlap |S^,(t)| = | J2 n (~ l)"0n(*)l> a l so shown in 
Fig. 4, becomes large and the CCR assumption breaks 
down and yields a qualitatitely wrong long time behaviour 
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with free acceleration as in the continuum. It should be 
pointed out that for a linear potential using a 27r-periodic 
funtion of ak for the kinetic energy, like e.g. the nearest 
neighbour hopping form (1 — cos ak)/ a 2 the CCR result 
for x(t) accidentally agrees with the exact solution. This 
is related to the fact mentioned earlier that the wrong 
use of the CCR yields the exact commutation realation 
[x, e ~ mak ] = ane~ mak for integer n. 

After this discussion of the linear potential we can re- 
turn to the harmonic potential and explain its properties 
in the regime where the CCR assumption breaks down. 
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ready rather small and the slow oscillation towards ana 
cannot be seen in the figure. Only in the very short 
time limit the CCR result for expectation value is useful. 
The fact that the deviations get more pronounced Jsv in- 
creasing no is opposite to what is widely believedtJ but 
is obvious from Fig.l. One obtains qualitatively similar 
results if the long range hopping is replaced by near- 
est neighbour hopping, i.e. k 2 /2 — * (1 — cos ak)/ a 2 (see 
dashed-dotted curve in Fig. 5), which corresponds to the 
Josephson Hamiltonian. The two regimes indicated in 
Fig. 1 correspond toJjhe Josephson and the Fock regime 
for this Hamiltoniariii In the latter regime assuming the 
CCR must certainly fail for the long time dynamics. 

In this paper we have shown that the fact that a canon- 
ical operator which is conjugate to the position operator 
on the lattice can be defined on a dense set of the Hilbert 
space is of very restricted use in practical applications. 
This is an indication that one has to be extremely care- 
ful in similar use of the particle number-phase CCR in 
fields like interactinjg_.electrons in one dimensionoQ and 
the Josephson effecirB. 

The author wants to thank V. Meden and W. Zwerger 
for stimulating discussions. 



FIG. 5. Average position of initially Gaussian wave pack- 
ets (b=0.2) centered at xq = — noa in a harmonic potential 
(c = 0.01) as a funtion of \fct. For no = 20 (dotted curve) 
the motion is harmonic as predicted by the CCR. For no = 30 
(dashed curve) many different excitation energies contribute 
and the motion is no longer harmonic. For no = 40 (full 
curve) the time scale is too short to show the slow oscillation 
towards noa and back. The dashed-dotted curve is for the 
case of nearest neighbour hopping and no = 20. 

The almost degeneracy of pairs of eigenvalues shown in 
the upper right corner of Fig. 1 can now easily be ex- 
plained. As shown in Fig. 3 (circles) the eigenstates 
for large enough quantum numbers are just even and 
odd states which (apart from a factor y2) almost look 
like the Wannier-Stark state which one obtains when the 
harmonic potential is linearized locally. The bonding- 
antibonding splitting decreases exponentially with the 
distance of the the right and left part of the eigenfunc- 
tions. This leads to a very slow component in the time 
evolution of the position operator shown in Fig. 5. For 
the same potential the time evolution strongly depends 
on the initial position xq — —noa of the wave paket. 
For sufficiently small values of no the expectation value 
(x(t)) is to a very good approximation given by a cosine 
function as predicted by assuming the CCR. For no = 30 
many different frequencies contribute and lead to clear 
deviations from harmonic motion (dashed curve). For 
no = 40 the initial state has strong overlap with eigen- 
states for which the bonding-antibonding splitting is al- 
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